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Abstract 

Several physical quantities of light hadrons are examined by a new holographic 
model of QCD, which is the modified version of the one proposed by Erlich et al. 
defined on AdSs. In our model, AdSs is deformed by a non-trivial bulk scalar, 
this is corresponding to adding the mass term of the adjoint fermions to the 4d 
SYM theory dual to the gravity on AdSs. We find that this deformation should 
be taken to be rather small, but its important effects are also seen. 
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1 Introduction 



Recently, the gravity/gauge correspondence has revived the expectation that QCD can 
be described by a string theory with an appropriate combination of D-branes. Some 
models are shown in the system of D p /D p+4: branes 0-0 • In these, by setting the 
lOd background with stacked D p branes, D p+4 branes are embedded as a probe to 
introduce flavor quarks. Then several physical quantities of QCD have been obtained 
with sufficient values. 

When these models are set as 5d theories, we should prepare two kinds of 5d ac- 
tions, Sbuik and S mesori . The former describes the bulk background or closed strings 
and the latter the hadrons or the system of open strings. They are corresponding 
to lOd action of gravity and the one of -D p +4 branes, respectively. The solution of 
SWik gives the gravity-background dual to the Yang-Mills theory. On the other hand, 
the classical solution of S meson provides parameters related to quarks, and the meson 
spectra are obtained from fluctuations around this solution in 5 mcson . However, there 
is no satisfactory holographic model, which describes QCD, at present. Instead, sev- 
eral authors have proposed phenomenological 5d-models which explain a wide range 
of hadron properties related to QCD EH HH H2J- In these models, the gravity- 
background is taken as AdSs and this space is cut off at an appropriate infrared (IR) 
point in order to realize the quark-confinement of the dual gauge theory. However, true 
gravity-background dual to QCD should be largely deformed from AdSs in the infrared 
region since the AdSs background is dual to M = 4 supersymmetric Yang-Mills (SYM) 
theory. Therefore, even if the infrared region is removed by the IR cutoff mentioned 
above, we expect that the deformation of the background configuration from AdSs is 
observed and it would affect the physical quantities. 

Our purpose is to make clear this point in the approaches mentioned above. This 
would be important to proceed these approaches to the next step. 

Here Smeson is set as the form used in 0. In this action, the mass and the vacuum 
expectation value (VEV) of bilinear fields of quarks are given through the classical 
solution of a tachyonic scalar field. As mentioned above, we deform AdSs by introducing 
the mass of adjoint fermions into the dual gauge theory. This is realized by adding a 
non-trivial scalar, whose conformal dimension (A) is three, to Sbuik- As a result, in 
the dual gauge theory, the supersymmetry is explicitly broken and the gauge coupling 
constant is running. So a new parameter is introduced in our model as the mass of the 
adjoint fermions, through which we examine its dynamical effects on various physical 
quantities related to light mesons. 

In the next section, we give our bulk-background configuration. In section 3, S meson 
and the configuration of the scalar of A = 3 are given, and physical quantities of light 
mesons are examined to see effects of our deformation. In the final section, summary 
is given. 

2 Bulk action and deformed SYM 

Up to now, many interesting 5d supergravity solutions have been studied as deformed 
SYM theories. Here we adopt a simple but nontrivial model given in [TH], which is 
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briefly reviewed. Consider the following 5d action with a scalar (0), 



S bulk = J d 4 xdz^ ^(R - 2A) - i(<90) 2 - , (1) 
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V( ( j>) = -%smtf{d ■-<{>} . i2) 



where [i = y — A/6. The parameters k 2 and A denote the five- dimensional gravitational 
constant and the cosmological constant, respectively. From (J2J), we can see that the 
mass of is M 2 = — 3yU 2 , then it corresponds to the conformal dimension three (A = 3) 
operator of TV = 4 SYM theory [T4] . 

We can solve equations of motion for metric and <p under the ansatz, = <p(z) and 

ds 2 = g M Ndx M dx N = A 2 {z) (ri^dx^dx" + dz 2 ) . (3) 

Here {M, N} and {/i, u} denote the 5d and 4d indices, respectively, and our metric 
convention is tjmn = ( — h + + +)• Setting as 

K 2 a 2 
A = ~ 

where a is a constant, we obtain the following solution 

M*)= .j.i . *(*) = as. (4) 
/i sinn VA^ 

From the second equation of (J3J), we can see that the parameter a represents the mass 
of the adjoint fermions of M = 4 SYM theory ^3] since <fi corresponds to the operator 
of A = 3 as mentioned above. In this sense, the supersymmetry of the gauge theory is 
broken for non-zero a, and the gauge coupling would be running. 



3 Meson spectra 



To discuss the meson properties, we start with the 5d action of the theory with the 
background and (jH). The fields in the bulk we consider here are the gauge fields, 
Lm and Rm, and a scalar field $ whose VEV is connected to chiral symmetry breaking. 
In this letter, we focus on Nf = 2 flavors. Then $ transforms as a (2l, 2r). The action 
is 



"Smeson — J d xdz \J ~Q Tr 



1 



4<? 



tMN i D r>MN\ 
2 \-"MNJ-> + tt-MN-ti 



\D M $\ 2 -M||$| 2 



, (5) 



where the covariant derivative is defined as Dm& = 9m$ + ^m$ — i$>Rm and g is the 
determinant of the metric. Lm = L a M r a , where r a are the Pauli matrices and similarly 
for other fields. g 5 is the 5d gauge coupling. We define <3> = S , e wr °' r0 and \v(z) = (S), 



where S corresponds to a real scalar and tc to a real pseudoscalar (S 
under L <-» R). They transform as 1 + 3 under SU(2)y. 



S and 7r 



-7T 
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3.1 Scalar field 



chiral symmetry breaking: Firstly let us study v(z) in the case of our background. 
Here and hereafter we take \i = 1. We put M| = —3 following 0110], which corresponds 
in the CFT side to an operator with the conformal dimension A = 3, through the 
general relation A = 2 + J A + M|. Solving the bulk equation of motion for S, 



(Ck-Ml)S 







where D5 denotes the five dimensional Laplacian, we obtain ^H] 
1 



v(z) 



-\smhr\A(z)/V\) + A(z)JA 2 (z) + A 



(6) 



(7) 



Here m q and c should be identified with the quark mass (explicit breaking of the 
chiral symmetry) and the chiral condensate (spontaneous breaking of chiral symmetry), 
respectively. Actually this is approximated as v(z) ~ m q z + cz 3 near z — 0, which 
reproduces the results of [TU| . We notice that v(z) includes the parameter A which is 
specific to our model as well as m q and c. In other words, the gauge theory considered 
here is characterized by these three parameters. 

d-meson: As for the singlet meson state (a), it is obtained as a solution for the 
fluctuation of S, S = v(z)/2 + a, with finite four dimensional mass m 2 defined as 
m 2 a. The equation for a is explicitly given in the same form with Eq. © as 
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A 2 



m 



+ d 2 + 3 



dzA 
A 



a 



Mia. 



It should be noticed that this equation is independent of v(z) and depend on M| and 
A through the warp factor A. 

Also in this case, the bound state spectrum is obtained by introducing an IR cut- 
off, z m , into the fifth dimension [31101, which corresponds to Aq CD . In this case, the 
boundary conditions are adopted such that a(z)\ Zo = d z a(z)\ Zm = 0, where z is the UV 
cutoff which is taken to zero after all. Then the mass of a depends on the parameters, 
A and z m since Mj is fixed. The numerical results are shown in the Fig^B), where the 
spectrum of vector fields are also shown in (A). In order to reproduce the experimental 
value m CT = (400 — 1200) MeV, A should be smaller than 0.35. The upper-bound for 
m a is shown by the dotted line in the Fig|TjB). 



3.2 Vector meson 

As for the gauge bosons, they are separated to the vector and the axial vector bosons 
Vm and Am, and are defined as L M = Vm + A M and Rm = Vm — A M , respectively. 

First of all, we discuss some properties of the vector mesons. The linearlized equa- 
tion of motion for is (employing V z — gauge) 

+ A- x d a M a \R = 0, (9) 
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where the mode expansion V^(x,z) = ^2 n ^u:H x )fn( z ) * s applied and is the four 
dimensional mass of the n-th excited vector meson* . The boundary condition for vector 



mesons is given by fX(z ) = d z f„ 



cVi 



0, similarly for the case of the sigma meson. 



Note that Eq.® is independent of v(z). Then the masses of vector mesons are given 
as the function of A and Z^yi £LS 111 the case of the a meson. 

It is convenient to solve Q numerically with the above boundary conditions. The 
first excited vector meson mass, m p , is obtained as the function of A and z m . By using 
the experimental value of m p , the IR cutoff z m is expressed as the function of A. As the 
result, the second and all highly excited vector meson masses depend on A. We find 
that the lowest and the next excited m p -meson are fitted by the parameters z m = 4 
and A = 0.8. 

The vector meson decay constant is computed through the second derivative of its 
own wave function according to jU] 



gl 







_ dz 2 


ZQ- 



gl 



12tt 2 



(10) 



Here N c = 3. However, its numerical value given for z m = 4 and A = 0.8 deviates from 
the experimental value. So we must change the value of A to the smaller side in order 
to obtain more reliable result. 




Fig. 1: (A) A dependences of the second excited vector meson m V2 (a) and the p meson 
decay constant F^ 2 (b) divided by their experimental values jTSj, which are used 
hereafter. (B) The curve and the dotted line denote m a (MeV) and its experimental 
upper bound, respectively. 



Our results are shown in Fig. Q where we find that a good fit would be obtained 
near A ~ 0.8 if we consider only the vector meson sector. Note here that in this sector, 
v(z) (or equivalently, the quark mass m q and/or the chiral condensate c) does not 
appear. This is contrasted with the results from, for instance, the QCD sum rules 
[OH El , where the vector meson masses depend on the chiral condensate. This is due 
to the fact that the spontaneous chiral symmetry breaking is introduced by hand in 
this model, which was done by the choice of the profile of the scalar field S. In this 
sense, the analysis can not be complete at this stage. 

* fn ( z ) i s the wave- functions for vector mesons with the normalization condition 1 = 
J z z -dzA(z)(fX(z)) 2 . 
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3.3 Axial-vector meson and 7r-meson 



The linearized equations of motion for the axial vector meson and the pion 7r are 
obtained from the following quadratic action 



S : 



axial 



dx dz 



A 



A 3 



{F a A Y - ^—{d<K a + A 



a\2 



where v(z) is given by (J7J). Decomposing A^ into the transverse and the longitudinal 
part, A^ = A^± + d^ip, one can obtain equations for these fields, 



[m\ + A~ l d z Ad z - glA 2 v 2 }A^ 
d z (Ad z v)-g 2 5 A 3 v 2 (7r + <p) = 0, 
m 2 d z <p + g 2 5 A 2 v 2 d z n = 0. 



0. 



(12) 
(13) 
(14) 



These equations are solved numerically under the boundary conditions, A fl ±(z ) = 
dzA^^Zm) = 0, ip(z ) = d z ip(z m ) = n(z ) = 0. The decay constants of the axial 
mesons and the pion are calculated from the wave functions as 



1 



\d 2 f* 




_ dz 2 


z . 



(n + 0), 



F 2 = - 



1 dj 



gl 



(15) 
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where A^ ± (x, z) = £ n A^x)f^(z) and f A (z) is normalized as 1 = f*™ dzA(z)(f£(z)) 2 , 
and / is the solution to Eq. (|12|) with m 2 A = 0, satisfying / = at z = z and df /dz = 
ciij z — Zfji • 

The masses and the decay constants of axial-vector mesons and pion depend on 
four parameters m q , c, z m and A through v(z) and A(z) in Eqs. (fT2^ - ([Ti|l . while those 
of vector mesons do only on z m and A. For the consistency between the vector and 
axial-vector meson sectors, here we take the same z m as that determined in the vector 
meson sector. The z m depends on A as l/z m = 0.323 — 0.09125A [GeV] for A < 0.8. 



m 2 [MeV 2 ] 




Fig. 2: The m q dependence of m 2 and m A1 . The solid line represents m 2 , and the 
dotted line corresponds to m 2 Al /100. The GOR relation is shown by closed circles. 



First, we show that the present model reproduces the Gell-Mann-Oakes-Renner 
(GOR) relation, m^F 2 = 2m q c, well satisfied in real QCD. Figure 121 shows the m q 
dependence of m 2 , where other parameters are fixed at A = 0.1, c = (0.3256 GeV) 3 
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and l/z m = 0.315 GeV. The solid line is a result of direct calculations of m 2 , and 
the closed circles are a result obtained from calculated F n through the GOR relation 
m 2 = 2m q c/F 2 . The two results agree with each other. Thus, the GOR relation is 
satisfied for the case of A = 0.1. This is also true for other A, as shown later. In 
this figure, when m q = 2.41 MeV, calculated m n and F n reproduce the corresponding 
experimental values simultaneously. For comparison, we also plot the m q dependence 
of m A1 by the dotted line. The m q dependence is quite weak. This means that the 
value of mA\ is determined not by m q but by the chiral condensate c. 

Figure El shows the c dependence of VfiA x and F Ai 2 , in which other parameters 
are fixed at A = 0.1, m q = 2.41 MeV and l/z m = 0.315 GeV. The two quantities 
in the axial-vector sector have a similar c dependence. They monotonously increase 
as c increases, and they almost reproduce the corresponding experimental values at 
c ~ 0.025 GeV 3 . Mass m Al is reduced to about 60 % of the experimental value in the 
limit of no chiral condensate. Thus, 40 % of the observed mass is generated by the 
finite chiral condensate. 




0.01 



0.02 



0.03 



c [GeV J 



1/2 

Fig. 3: The c dependence of F A and m Al . These quantities are devided by the 
corresponding experimental values. 



We focus our analysis on the A dependence of F Al , m^, F w . Parameters m q and 
c are fixed as follows. First, c is assumed to be determined from m q through GOR 
relation, rh 2 F 2 = 2m q c, where fh n and F n denote experimental values of m n and F n . 
Second, m q is fixed so as to reproduce the observed pion mass. The resultant m q 
depends on A as m q = 2.26 + 1.125A [MeV] for A < 0.4. 

nv 1 /2 

Figure EfA) shows A dependence of predicted values of F A ) mA 1 -,F 7T . As for all A 
up to 0.4, calculated F n agrees with the experimental value with ~ 1% error. In the 
present analysis, we assumed the calculated F^ to reproduce the experimental value 
when we used the GOR relation as a relation between c and m q . The good agreement 
shows that our assumption is consistent and then the GOR relation is well satisfied for 

1 /2 

all A at least up to 0.4. As for F A and m Al , the agreement of the theoretical results 
with the observed ones becomes better as A decreases. One then see that the A = 
case, namely the AdSs case, yields a best fit if we consider only the axial-vector meson 
sector. 

Figure EfB) summarizes all predictions of our model in both the vector and axial- 
vector meson sectors. As for three of five quantities, F A± 2 ,m Al , F^ 2 , the agreement of 
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Fig. 4: Predictions of the model in the axial-vector sector (A) and in both the vector 

1 /2 

and axial- vector meson sectors (B). As for (A), the solid lines, a, b, c, denote FJ , 

1 /2 

m^, F n , respectively. As for (B), the solid lines, d, e, g, h, represent Fj[ , m^, F n , 
F^ 2 , respectively, and the dotted line, f, corresponds to m V2 . All quantities are divided 
by the corresponding experimental values. 



the predictions with the corresponding observations is best at A = 0, but as for one of 
the five, i.e. for my 2 the best fit is obtained at A ~ 0.8. Hence, we can conclude that 
totally the best fit is realized at small A. This result is welcomed to obtain a realistic 
mass of a meson shown in Fig^ 



4 Summary 

We present a new holographic model of QCD in which the gravity-background is de- 
formed from pure AdSs and the extra dimension is cut off at an appropriate infrared 
point. This is a natural extension of the holographic models of Refs. EH EH E2] 
in which the AdSs background is taken and the extra dimension is cut off in the same 
manner. These holographic models based on AdSs well reproduced observed quanti- 
ties on light mesons. Nevertheless, there is no strong reason why AdSs is taken as a 
background of the gravity dual to QCD, since QCD has a running coupling constant 
and then the gravity-background should be modified from AdSs. Our new model is 
proposed to answer this question. 

In our model, parameter A represents the magnitude of the deformation of AdSs 
background. For any A, our model reproduces the GOR relation which is well satisfied 
in real QCD. As for light mesons except the excited p meson, a best fit of our predictions 
to the corresponding observed quantities is obtained at A = 0, that is, in the AdSs limit. 
Meanwhile, the excited p meson mass is reproduced at A = 0.8. Hence, a best fit to all 
of them is realized at small A. The parameter A is related to the five-dimensional Planck 
mass M 5 and the scale of the supersymmetry breaking Asusy as A = A| USY /(3Mf). 
The smallness of A would imply M 5 ^> Asusy i n the gravity dual to QCD. 

It is well known that the chiral condensate c plays an important role on all quantities 
of light mesons. In this sense, the quantities should depend on c. Our model has the 
property for the axial-vector meson sector, but not for the vector meson sector. Thus, 
the present model is insufficient for the vector meson sector. The holographic models 
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of Refs. [3 El El have essentially the same problem. This is an important problem 
to be solved in future. 



Acknowledgments 

M. T. would like to thank T. Hatsuda and D. T. Son for valuable information. This 
work has been supported in part by the Grants- in- Aid for Scientific Research (13135223) 
of the Ministry of Education, Science, Sports, and Culture of Japan. N.M. was sup- 
ported by RIKEN (No.Al2-61014). 

References 

[1] A. Karch and E. Katz, JHEP 0206, 043(2003) |hep-th/0205236 . 

[2] M. Kruczenski, D. Mateos, R.C. Myers and D.J. Winters, JHEP 0307, 049(2003) 
|hep-th/0304032| . 

[3] M. Kr uczenski, D. Mateos, R . C. Myers and D. J. Winters, JHEP 0405, 041 
(2004) |arXiv:hep-th/0311270 . 

[4] J. Babington, J. Erdmenger, N. J. Evans, Z. Guralnik and I. Kirsch, Phys. Rev. 
D 69, 066007 (2004) |arXiv:hep-th/0306018| . 

[5] N . J. Evans and J. P. Shock, Phys. Rev. D 70, 046002 (2004) 
|arXiv:hep-th70403279|. 



[6] C. Nunez, A. Paredes and A.V. Ramallo, JHEP 0312, 024(2003) |hep-th/0 311201|. 

[7] T. Sakai and S Sugimoto, Prog.Theor.Phys.ll3(2005)843-882, lhep-th/041214l"| ; 
|hep-th/0507073| 

[8] K. Ghoroku and M. Yahiro, Phys. Lett. B 604, 235 (2004) | arXiv:hep-th/0 408040|. 

[9] J. Erlich, E. Katz, D. T. Son and M. A. Stephanov, [hep-ph/ 05 011281 

[10] L. Da Rold and A. Pomarol, Nucl. Phys. B 721, 79 (2005) |hep-ph/ 05 012181 ; 
|hep-ph/0510268| 

[11] G. F. de Teramond a nd S. J. Brodsky, Phys. Rev. Lett. 94, 201601 (2005) 
|arXiv:hep-th/0501022| . 

[12] J. Hirn and V. Sanz, |hep-ph/0507049 

[13] K. Ghoroku, M . Tachibana and N. Uekusa, Phys. Rev. D 68, 125002 (2003) 
lhep-th/0304051"] . 

[14] L . Girardello , M. Petrini, M. Porrati an d A. Zaffaroni, J HEP 9812 (1998) 022 
|hep-th/9810126| ; JHEP 9905 (1999) 026 |hep-th/9903026| . 

[15] K. Ghoroku and M. Yahiro, Class. Quant. Grav. 20, 3717 (2003) |hep-th/021lTT2 . 

8 



[16] N. V. Krasnikov, A. A. Pivovarov and N. N. Tavkhelidze, Z. Phys. C19, 301 (1983). 

[17] L. J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rept. 127, 1 (1985). 

[18] S. Eidelman et al. [Particle Data Group Collaboration], Phys. Lett. B592, 1 
(2004). 



9 



